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Abstract. We investigate an action that includes simultaneously original and 
dual gravitational fields (in the first order formalism), where the dual fields 
are completely determined in terms of the original fields through axial gauge 
conditions and partial (non-covariant) duality constraints. We introduce two kinds 
of matter, one that couples to the original metric, and dual matter that couples 
to the dual metric. The linear response of both metrics to the corresponding 
stress energy tensors coincides with Einstein's equations. In the presence of 
nonvanishing standard and dual cosmological constants a stable solution with a 
time independent dual scale factor exists that could possibly solve the cosmological 
constant problem, provided our world is identified with the dual sector of the 
model. 
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1. Introduction 

The known concept of Hodge duality for form fields (antisymmetric tensor 
fields), i.e. the permutation of free equations of motion and Bianchi identities, can 
be generalized to fields of mixed symmetry [l]-[3] . Here again, linear combinations 
of free equations of motion and Bianchi identities imply Bianchi identities and free 
equations of motion for the dual fields or their field strengths. In particular, this 
generalized Hodge duality can be applied to linearized gravity in d = 4 or higher 
dimensions [T]-|ll|. 

One motivation for the introduction of a dual graviton are attempts to realize 
hidden symmetries in d = 11 supergravity/M-theory [T2]. However, only in d — 4 
the dual of a graviton is again a symmetric two component tensor field pQ- [9], that 
could be used to construct a dual gravitational theory. 

Any attempt in this direction has to face no-go theorems on local Lorentz 
invariant dual models of gravity [13]-[14j. The way out considered here is not to 
insist on a duality symmetry, and - most of all - to break Lorentz symmetry twice: 
Once through a non-covariant (axial) gauge fixing for the dual gauge fields, and 
through non-covariant duality constraints on the field strengths. 

Non-covariant gauges have the particular property that, in the context of Yang- 
Mills theories, the gauge field can be reconstructed from the field strength [151 IS] • 
This feature has been used in [15] in order to construct dual models of Yang-Mills 
theories that are not, however, duality symmetric. 

Here we study a similar program in the context of 4c? gravity. In order to 
formulate gravity as closely as possible to Yang-Mills theories, we employ the 
first order Cartan formalism where vierbeins and (spin) connections are treated 
as independent fields, and the original local symmetries are diffeomorphisms and 
an "internal" 0(3, 1) Lorentz symmetry (the gauge fields of the latter are the spin 
connections). 

Let us first outline the rough idea behind the present approach: the starting 
point is a standard gravitational action involving vierbeins e a ^ and connections u ab fl 
(here and in the following Latin indices indicate representations under the internal 
Lorentz group, whereas Greek indices are space-time indices). 

The Riemann-Cartan tensor R ab ^ u is the field strength of the connections, and 
its dual R ab ^ u can be constructed by the contraction of a pair of (antisymmetric) 
indices with two of the four indices of the 4d completely antisymmetric epsilon 
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tensor. 

Once R n V is given through a duality relation in terms of R ab ^ ul one can 
ask under which circumstances R ab lll/ can be written as a field strength of a dual 
connection uj ab ,,. 

First, for each fixed choice of [a b], R ^ u contains 6 antisymmetric combinations 
of indices [fi v\. Deducting the gauge degrees of freedom from 0J ab fl , we have however 
only three independent connections uj ab fl (for given [a b}) at our disposal. This leaves 
us with two possibilities: (1) R ab ^ v has to satisfy certain constraints; (2) we only 
use 3 out of the 6 antisymmetric combinations of indices v\ in R ab fiu in order to 
construct u) ab „. 

The constraints (1) would be the (three independent) Bianchi identities. In 
fact, in a weak field expansion the validity of the Bianchi identities for R ab ^ u can 
be deduced from the ones for R^^ under the condition that the Ricci tensor R^ v 
vanishes Beyond a weak field expansion (where the covariant derivatives 

acting on R al ' would involve the not yet known connection u) ab ll , that differ in any 
case from uJ ab fJ ), or for no n- vanishing R^ u , the Bianchi identities for R ab „ u can no 
longer be employed. We are left with the possibility (2) above: we construct u; a6 
from the 3 combinations of indices [fi v\ that are left over after a contraction of 
R a> with a fixed vector n v '. Note that we still assume that all 6 combinations of 
indices \p v\ of R ab ^ v are identified with the field strength of the connection cD ab M ; 
we do not assume, however, that the 3 additional combinations of indices [p, v\ of 
R al ' are determined through a duality relation in terms of -R 06 ^: 



First, for a fixed vector nn§|, we define only 3 combinations of indices [/i v\ of 
R ab ^ in terms of R ab ^. 

n"R a \ u = \ e ab cd n»R a \ u . (1.1) 

Next we introduce a dual connection u> ab fl . On the dual connections we impose 
the axial gauge condition 

u ab ^ n» = . (1.2) 

Then we interpret R b ^ y as a field strength for u , which implies - in the gauge 
f OD - the relation 

rfdp u ab u = n»R ab ^ . (1.3) 

§ In section 2 we actually employ a Lorentz vector n a , but the essential idea is more easily 
formulated with n^. 
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Once n^R ab ^ u is a known expression in terms of R ab ^ u , the enforcement of both 
conditions (11. 2p and (11.31) determines the dual connections completely (up to 
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boundary conditions on a hypersurface perpendicular to n^) in terms of R 

A similar reasonning can be employed for the construction of a dual vierbein 
e° , out of which a dual metric g^ u = e a ^e a u can be constructed: For given u) ah ^ the 
dual vierbein could be constructed from the vanishing of the torsion tensor 

T;, = %e^ ]+ ^e% . (1.4) 

Again - for each fixed value of a - T a ^ u contains 6 antisymmetric combinations of 
indices [/x u], but (modulo gauge degrees of freedom) we have only 3 degrees of 
freedom in e a M at our disposal. Again we can restrict ourselves to the vanishing of 
T a n u ', which allows to construct the vierbeins e a (in the axial gauge), but the 3 
extra components of T a ^ u (for each a) will vanish no longer in general. (One can 
prove that they would vanish in a weak field expansion, but our aim is to go beyond 
it.) 

Hence, we impose a gauge condition on e a ^ (i.e. a choice of the coordinate 
system) such that 

Z\ n» = 8\ n» (1.5) 

where 5°"^ is the Kronecker symbol. A coordinate system where ( jl.5p holds can 
always be constructed locally. 

Contracting (jl.4p over and using both eqs. (11.2j) and (11.51) one obtains 

n"0 M ?„ = u\ v S b X . (1.6) 

Similar to the case of the connection, eqs. (11.51) and (11.61) determine the vierbeins 
e 0, completely up to boundary conditions, now in terms of the connections a;" 6 
constructed before. Again, eqs. (11.51) and (11.61) do not imply the vanishing of all 
components of T a M „ (again only half of them), i.e. the dual theory is not entirely 
torsionless. 

On the other hand we have now achieved our goal: the introduction and 
definition of dual gravitational fields, including a dual metric g^, beyond a weak 
field expansion. Moreover, all the previous conditions (jl.ip - (11.31) . (11.51) and (II. 6p 
can be cast into the form of an action involving Lagrange multipliers. 

The possibility to construct an action for dual, and generally interacting, fields 
is highly non-trivial: First, attempts to impose (covariant) duality relations for 
all components of field strengths using Lagrange multipliers leads generically to 



equations of motions for these Lagrange multipliers that have non-trivial solutions 
corresponding to new degrees of freedom. Attempts to eliminate all new degrees of 
freedom through new multipliers result in the need for an infinite tower of Lagrange 
multipliers [17J. The corresponding problem can be cured quite easily here, with 
just one extra multiplier required (see section 2). 

Second, the complete permutation of non-linear equations of motion and 
Bianchi identities (between original and dual fields) resulting from an action is 
generally impossible, since an action allows for the addition of sources for fields 
that modify the equations of motion, but leave Bianchi identities unchanged. Here 
duality relations are satisfied by only some of (half of) the field strengths - that 
satisfy the second Bianchi identities by construction - and, furthermore, the first 
order formalism allows for torsion that allows to violate the first Bianchi (cyclic) 
identities of the Riemann tensor. 

In fact, the possibility to construct an action involving simultaneously the 
original and the dual gravitational fields as in the following section 2 allows in 
particular to couple the dual metric, in a standard fashion, to matter and to study 
its reactions to a stress energy tensor. We will carry out this analysis in the weak field 
limit in chapter 3, with the result that this reaction is the same as in standard general 
relativity. This implies also that vacuum solutions of standard general relativity 
remain valid also for the dual metric in the weak field limit, but we have to expect 
modifications of the solutions of the combined set of equations of motion (for both 
the original and the dual fields) beyond lowest order. In particular any attempt 
to integrate out completely the original gravitational fields will result in non-local 
interactions for the dual fields, hence their corresponding effective theory differs 
definitely from standard general relativity. 

The present approach is manifestly asymmetric between the original and the 
dual gravitational fields; it seems generally impossible, however, to implement a 
duality symmetry into interacting gravitational theories [T3~t [T4"]. 

In section 4 we study cosmological solutions of the equations of motion derived 
in section 2. We consider simultaneous cosmological constants in the standard and 
the dual sector and obtain a stable solution with a static dual scale factor. This 
scenario could provide a solution of the cosmological constant problem, provided 
our world is identified with the dual sector of the model. 

Finally, a summary and outlook is given in section 5. 
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2. The action and its variation 



The basic idea for the construction of an action involving simultaneously 
original and dual gravitational fields has been outlined in the introduction: Dual 
gravitational fields (the connection and the vierbein) are introduced together with 
constraints that determine them completely in terms of the original fields, up to 
boundary conditions. This procedure guarantees the absence of new degrees of 
freedom propagating in 4d space time, but the price to pay is the explicit breaking of 
both internal (Lorentz) and space-time (local coordinate reparametrizations) gauge 
symmetries. (The fixation of space-time gauge symmetries could, in principle, be 
omitted, but then the absence of new degrees of freedom is less obvious.) 

In this section we present an action, and discuss the various steps of symmetry 
breaking. In order to keep track of local coordinate reparametrizations it turns out 
to be convenient to proceed slightly differently than outlined in the introduction 
(although the final result is the same): Instead of introducing a fixed constant 
vector n M , that would not be invariant under local coordinate reparametrizations, 
we introduce a fixed constant Lorentz vector (but space-time scalar) n a that is not 
invariant under local or global Lorentz transformations. 

We start with the standard Lagrangian for general relativity in the first order 
Cartan formalism, where the independent fields are the connection u ab fl and the 
vierbein e a ^ (or its inverse ej 1 ). Lorentz indices a, b, ... are raised and lowered with 
the flat Lorentz metric 

?U = diag (1,-1,-1,-1) . (2.1) 
The Riemann-Cartan curvature tensor is the field strength of the connection: 

R a \ v = dp u a \ - d v co a \ + W % u/*„ - oo a cu oo c \ 

= 2(^^ ]+ < [M ^) (2.2) 

The inverse vierbeins allow to construct a space-time scalar version of the Riemann- 
Cartan curvature tensor, 

R ab cd = e/ e/ R a \ u . (2.3) 
Then the standard Einstein Lagrangian is 

C E = ^det(e) R ab ab (e,cu) . (2.4) 
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(Our convention here and below is that all terms in the Lagrangian transform as 
densities under local coordinate reparametrizations, such that the final action is 
simply / d A xC). 

Next we introduce a new "dual" connection oj ab ,,, and a new "dual" vierbein e a „ 
(or its inverse e a ). At the first level concerning tu ab ^, we have to add three Lagrange 
multiplier fields: 

i) K will serve to impose the axial gauge condition on uj ab ^ 

ii) LP^Jj will serve to impose that the dual of the Riemann tensor (12.31) . 
contracted with n d over its last index, is equal to the field strength R ab cd of the 
dual connection uj ab ^ (defined in analogy to (12. 2p ). again contracted with n d . 

A few comments are in order here: In the case of standard duality between 
(abelian) antisymmetric tensor fields, the dual field strength is obtained by 
contracting an epsilon tensor (with space-time indices) over the space-time indices 
of the original field strength, which would correspond to the lower pair of indices 
of the Riemann tensors here. Below we need, however, one of the lower indices - 
contracted with n d - in order to employ a relation similar to eq. (11.31) . Therefore 
we employ the upper pair of indices in order to define a dual Riemann tensor. 
In the absence of torsion, this would make no difference, since then the Riemann 
tensor would be symmetric with respect to an exchange of these two pairs of indices. 
Below we will find, however, that torsion is generally non-vanishing here, i.e. not all 
components of the connections 0J ab ^ and are related in the standard fashion to 
the corresponding vierbeins. Hence the present definition of a dual Riemann tensor 
in the sense of duality with respect to the internal 0(3, 1) Lorentz symmetry differs 
somewhat from standard S'-duality. 

Finally we found it convenient to equate the Riemann tensors multiplied with 
the determinants of the corresponding vierbeins, which simplifies some of the 

(1) c 

equations of motion below and allows L ab to transform as a space-time scalar. 

iii) It turns out that the components of LP~\£ with the index c in the direction 
of n c , i.e. n c La£, are not constraint by the action. In order to eliminate these 
unwanted degrees of freedom we employ another Lagrange multiplier N^ ab . No 
further Lagrange multipliers are needed in order to cope with unwanted degrees of 
freedom of N^ ab . 

Hence at this first level, which treats the constraints on the dual connection 
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w ab ^, the following three terms are added to the Lagrangian: 



+ L {1) J n d ^det(e)e ab ef R e Ue^) - det(e)R ab cd (e,ti)} 

+ N^ ab n c L (1) a6 c . (2.5) 

(Here K^ l > and A^ 1 ) transform as densities like det(e) under coordinate reparamet- 
rizations) . 

For a given dual vierbein e a , all components of tD ab M are now fixed by the 
constraints following from the variation of (12.5ft with respect to K ^ and L^\^. 

At the next level, we introduce constraints on the dual vierbein e a M (or, for 
convenience, on its inverse e a M ). Again we will not follow exactly the procedure 
outlined in the introduction (although the final result will be the same), since we 
employ the constant Lorentz vector n a rather than a space-time vector n M . 

Otherwise, the roles of the three Lagrange multipliers introduced at this second 
level ressemble to the ones in (12.51) : 

i) K( 2 ' v serves to impose an axial gauge condition on e; 

ii) L/( 2 ' a b serves to impose the vanishing of the torsion tensor T a bc contracted 
with n c , with 

T\ c = X T%, (2.6) 

and T\ v as in (PP . 

iii) since the components L^ 1 \ h n\ ) are not determined by varying the action, 

N {2)a 

serves to eliminate these degrees of freedom. 
Hence at the second level, the following three terms are added to the Lagrangian: 

£(2) = K {2) u n b Q y + ^ b + N (2)a ^JQ b^ ^ 

(Here and transform as densities under coordinate reparametrizations) . 

Finally we can couple both kinds of gravitational fields to matter, after defining 
the metrics 

9\iv = ^p,a &v (2.8a) 

9»v = e u a . (2.8b) 

For instance, one can introduce matter that couples to g^, and another kind 
of matter that couples only to g^ u ; hence two kinds of matter Lagrangians 

-C M (g)+£ M (g) . (2.9) 



(The relative minus sign in (12. 9p serves just to reproduce the standard signs in the 
corresponding Einstein equations, see section 3.) 

Next we will derive the constraints imposed by the various Lagrange multipliers. 

Starting with £W in (12.51) . its variation with respect to K^ ab imposes a gauge 
condition of the axial type on u) ab tl : 

n c ^uj a \ = n»uj a \ = 0. (2.10) 
Below we will find that for constant n a , the vector 

n^ = n a ef (2.11) 
is constant as well. 

The variation of with respect to L^ ab c includes, at first sight, a term 
proportional to N^ ab . However, contracting the variation with n c and using the 
antisymmetry of the curly bracket in (I2.5P in [cd], one obtains 

N m*b = _ (2.12) 

Furthermore we can use that, with the definition analogously to (12.21) and (12. 3p of 
R ab cd and fl2HDD, 

R ab cd n d = -V r?d v (2.13) 
Hence one finally obtains the constraint from the variation w.r.t. L ab c 

det(e) if n"d v = -\ det(e)e ab ef R ef cd (e, w)n" . (2.14) 
The variation with respect to N^' ab trivially implies 

L {1) ab c n c = 0. (2.15) 

Now we turn to in ( 12. 71) . with the torsion tensor T a bc = T" 6c ] defined in (12. 6\j 
and (11.41) . The variation with respect to K^ 1 ' gives 

d v (n b e 6 M ) = d v n» = . (2.16) 
Contracting the variation with respect to L^ 2 a b with n b gives 

iV (2)a = . (2.17) 

Using (I2.17P and contracting the variation w.r.t. L^ 2 \ b with e b 1 gives 

n" (d, e\ - d v ~e\ + D%, e\, - u\ v e 6 J = . (2.18) 

Using (I2.16P and (12.111) one finds that the first term in (I2.18P vanishes, whereas 
(I2.10p implies the vanishing of the last term. Hence (12.181) collapses to 

^d v e\ = n b Uj\. (2.19) 
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Finally the variation with respect to N^ a trivially implies 

H 2 l c n, = . 



(2.20) 



Let us study the consequences of the constraints (I2.16P and (12.191) on the dual 
metric g^ v in (12.8bj) . Contracting (12.191) with n a gives 

n v d v (n a e\) = . (2.21) 
From n^g^ = n a e a ^ one finds that (12.211) is equivalent to 

n v d v (n" g^) = . (2.22) 
Axial gauges of the form (I2.22j) - with constant n M as in (I2.16P - have been used 



in [TS] in order to construct graviton propagators and to analyze one loop diagrams 
in quantized gravity. They can constrain the global properties of space-time; the 
Schwarzschild solution in such a gauge - with n M time-like, goo = — 1, goi = 0, i.e. 
in a comoving frame (or Novikov coordinates) - is known, however, and given e.g. 
in [TJ]. In cosmology, on the other hand, this gauge is standard. 

Next we give the equations of motion that follow from the variations of the 
complete action with respect to the vierbeins and the connections. 

First, from the variation with respect to e M , contracted with ej 1 and after 
division by det(e) one obtains 



1 



£ a r)bc , ryab 
" 2 °g be i* 11 gb 



+ 



[c d] mn 
i ef 



--S a n R ef cd + 6 a r R 



yd 



rp a 

1 Mg 



where 



T 



5 



Mg 



det(e) 8ef 



C 



M 



(2.23) 



(2.24) 



The variations with respect to u are best expressed in terms of the two tensors 



and 



TP 

^ab 



q v 



ab] 



det(e) e b v 



S, 



ab] 



[/us] 



det(e) L%^s ef ab e-ef 



Then one obtains 
1 

K 

For S, » v 



u f ^ab 



2E„ 



CO. 



b] v 



9 C P u 



b] v 



(2.25) 
(2.26) 
. (2.27) 



ab — 0, eq. ( 12.271) would determine lo in terms of e a M , as usual, through the 
vanishing of the covariant derivative of e° , or the vanishing of torsion. Generically, 
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however, S ab tiV will not vanish (since L^j will not vanish), hence torsion is 
generically non- vanishing for configurations u) ab fl that solve eq. (12.271) . 

Concerning the variation with respect to e a M , it is convenient to consider the 
combination 

V < 2 ' 28 > 

that is independent of the Lagrange multiplier field v in (EZJ). (The remaining 
components just determine K^ 2 ' u , that appears nowhere else.) Then one finds 




(2.29) 



with 



~ m S ~ 

Tm * = de% Mf Cm ■ (2 - 30) 

Finally the variation with respect to u) gives, again after an elimination of 
K^ h through a contraction with n d ej 1 , 

n d d v (det(e) L^ ab c - \ L % c n b] = . (2.31) 

Some consequences of these equations of motions will be studied in the next sections. 

Before concluding this section we have to make some comments on the Lorentz 
symmetry breaking induced by the constant Lorentz vector n a . Apart from the 
terms that impose axial gauge conditions (the first terms in ( 12.51) and ( 12. 7ft ). n a 
appears in various other terms in (I2.5P and (12.71) . In the next section we study the 
equations of motion in a weak field expansion (around a Minkowski vacuum, and 
considering all Lagrange multipliers as weak fields). At this level we find no explicit 
Lorentz symmetry breaking in the equations of motion. 

Gravitational self-interactions, that start to play a role in higher order in a 
weak field expansion, will generally not respect Lorentz invariance, however. In 
the worst case this could induce violations of unitarity and/or strong gravitational 
self-couplings at inacceptable length scales. 

It may then be advisable to promote n a to a field n a (x), and to replace 
explicit Lorentz symmetry breaking by spontaneous Lorentz symmetry breaking via 
a potential X(n a n a ± l) 2 in the Lagrangian for n a (x). Apart from a massive radial 
mode, this scenario will add three Nambu-Goldstone (NG) modes to the model. 
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The possible fate of such NG modes has recently been reviewed in [20J, and 
depends heavily on the Lagrangian (e.g. additional R 2 terms, that generate a 
propagating spin connection): the NG modes could remain massless, but could 
also become massive degrees of freedom. Some phenomenological consequences of 
such NG modes have been studied in [21] , but even the quadratic part of an effective 
Lagrangian would have been to be constructed here. 

Also, unless the Lorentz gauge symmetry is restored this way, the axial gauge 
condition (2.5) imposed on u> ab are not just gauge conditions that can, in principle, 
be chosen at will (in contrast to the gauge conditions (2.7) on e*): observables will 
depend on the present choice of the gauge condition imposed on uj ab ^ - however, 
only to higher order in a weak field expansion, see the next chapter. 

3. The linearized gravitational field equations 

The aim of this section is to study the linearized response of the gravitational 
fields to matter sources, notably to a dual stress energy tensor T Ma b (g). 

The point is that, as long as T Ma b (g) does not depend on g and vice- versa (what 
we will assume in the following), standard matter described by T Ma b (g) "sees" a 
space-time geometry described by g^, whereas dual matter described by T Ma b (g) 
"sees" a space-time geometry described by g^ that will generally be quite different! 
An interesting question is then whether g^ v reacts to T M b in a way that resembles the 
standard Einstein equations; if this is the case, our world could possibly correspond 
to the dual matter of the model. 

An answer to this question is not quite trivial even to lowest order, where 
(counting Lagrange multipliers as weak fields) terms ~ L ■ R in eqs. ( 12. 23ft and 
(I2.29P can be neglected: whereas the effect of T M b on g^ u is obvious from (12.23P 
(and agrees with Einstein's gravity), the effect of T M b is 

i) a non- vanishing value of L j, c from (12.291) . which induces 

ii) a non- vanishing value of L^\ h c from (12.311) . which induces 

iii) non- vanishing torsion for u ab fl from (12.271) . hence non- vanishing components 
of R ah ^ v1 which induce 

iv) a non- vanishing u ab ^ from (12. 14ft . that finally generates 

v) a non- vanishing e a fl from ( 12.191) . 

Subsequently we will carry out these steps explicitly. First it is convenient, 
however, to adopt a convention concerning the direction of n a , that we assume to 
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be time-like: 

n a = (l,0,0,0) (3.1) 

and the first component will be denoted by 0, i.e. n° = 1. All space and Lorentz 
coordinates perpendicular to n M ~ n a (to lowest order) will be denoted by latin 
letters i,j, k, . . . from the middle of the alphabet. 

Next we have to comment the fact that eq. (12.291) contains a projector such 
that the components T M0 d do not appear. Such a situation is actually familiar from 
Yang-Mills theories in axial gauges: After imposing the axial gauge condition, the 
corresponding components of the currents (here: the stress energy tensor) decouple 
from the gauge field. They contribute nevertheless to the dynamics of the theory 
through the equations associated with (covariant) current conservation, that can be 
derived via Noether's theorem. Likewise, we have to use the covariant conservation 
of T Ma b (and its symmetry) here, which gives to lowest order 

do f Ma ° + d % f M J = (3.2) 

for all a. Once T is assumed to be symmetric, eq. (I3.2p determines all components 
of T Ma b in terms of T Mi J , hence the absence of T M0 d in ( 12. 29ft constitutes no longer 
a paradox. 



Now we turn to eq. fl 2 . 2 9 [) which reads to lowest order (where e d ^ ~ 5 
e u ~ 8" , and only b — j contributes) 



d L {2 Y = 2T M I . (3.3) 



Mj 

The index combination (ab) = (jO) of eq. (12.311) gives, with (13.31) for 

1 
2 



do d L ( V = t T M / , (3.4) 



whereas the combination (ab) = (jk) implies the vanishing of L^- k % up to terms 
linear in t. Eq. (13 .4p implies a non- vanishing value for in (12.261) . which gives 

C fcO _ JO r(l) k _ 10 1 fp k l o c\ 

°ij — fc ij ^ 10 ~ b ij 2d d 11 
where we allowed ourselves to represent integrals with respect to x° = t by an inverse 
derivative 1 /do- Next we have to determine u ab fl from (I2.27p . and it is convenient to 
decompose u; a £, M into 

uJabf, = Q ablx (e) + (3.6) 



where 



^W(e) = - (d b (e ail + e Ma ) - d a (e bfl + e^ b ) + d^ (e ba - e ab )) (3.7) 
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and a>J 6 represents torsion. If one would replace u by Q in (|2.27p . its first line would 
vanish identically. Hence (I2.27P determines u T in terms of S^" ', and Q (or e) has to 
be determined elsewhere (by eq. fl 2 . 2 3 [) ) . Exploiting the various index combinations 
of eq. (12.271) one obtains after some calculation 

l fc ik 1 Mjn~rt. ij ±Mkn~rt k j ±Min) j 



T 
U ijk 



2d 



i0j ~ ~2d~o ij °™ ' 

Next we wish to construct R ab cd (u) from R ab cd {oj) via (I2.14p or, better, directly 
from the constraint from (12.51) : 

R a \ d = \ e ab ef R ef cd . (3.9) 

At first sight we have a problem, since R e ^ cd depends on uj = Q + u T , and Q is 
not yet known. However, it turns out that the linearized dual Ricci tensor does not 
depend on Q, since e ab e fR e ^ ad (Q) = identically. This allows us to construct 

R\ = R ab cb = \e ab ef R eS cb ^ T ) (3.10) 

with (13. 8p for u T , and using (13.21) in order to rewrite spacial derivatives. The result 
is 

R\ = k(t m \- 1 -8\T m c ^ (3.11) 

which seems to coincide with the standard Einstein equations. 

However, R a is the Riemann-Cartan tensor defined in terms of u> , and 
coincides with the Riemann tensor R ab ^ v {e) only if the dual connection u) 0- ^ is 
torsionless, i.e. if 



= fi«» M (e) (3.12) 

with Q(e) as in (13. 7p . In order to study this question, we have to construct all 
components of uj ab /I from eq. (12.141) . and subsequently e a M from (I2.19p . The result is 
that (I3.12p holds indeed, provided 

i) the original Ricci tensor R a b (uj = fi + uj t ) satisfies 

R°i = o , 

<9 i^ = 0; (3.13) 
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ii) the off-diagonal components e° i of the dual vierbein, that are required to be 
t-independent from eq. (I2.2ip . satisfy 

<%e°,] = 0. (3.14) 

(Actually, we could have imposed n a e a fjL = n a 5" , i.e. e°j = 0, from the beginning). 

Equations (13.131) follow from the not yet considered eq. (12.231) if T Mg a = 0, but 
they also allow for a cosmological constant T Mg a = S g a A in the "standard" matter 
Lagrangian. 

Note finally that for a vanishing dual stress energy tensor, the equations (I3.8P 
above imply vanishing torsion for the standard spin connection, and the linearized 
eqs. ( I2.23P can be interpreted as standard torsionless Einstein equations. 

To summarize, we have learned about two important features of the present 
model in this chapter: 

i) the possibility to reproduce the linearized Einstein equations including matter 
for the dual gravitational fields, in spite of the Lorentz-non-covariant action, 

ii) the crucial role played by torsion within the present first order formalism: 
"standard" torsion is induced by "dual" matter (and vice versa); this allows to 
generalize the known correspondence between vacuum equations of motion and 
Bianchi identities pQ-[7] to equations of motion with sources, a possibility already 
advocated (in a cosmological context) in [22J. 

4. Cosmological solutions 

In this section we return to the full nonlinear equations of motion of section 2 
and study cosmological solutions. The aim is to check under which circumstances 
the standard Freedman- Robert son- Walker (FRW) equations are reproduced, but 
also to see whether the model could provide a hint for a solution of the cosmological 
constant problem (CCP). A solution to the CCP would correspond to a (stable) 
solution of the equations of motion where A and/or A are non- vanishing, but the 
standard and/or dual metric remains nearly time-independent (does not explode 
exponentially with cosmic time). Indeed we will find such solutions for the dual 
metric, with non-vanishing A and A, below. 

First we have to make a general ansatz for all fields and Lagrange multipliers of 
the model, that is consistent with a homogeneous and isotropic universe (i.e. depend 
on x° = t only). For the Lorentz vector n a we will make the same choice as in eq. 
(13.11) . Note that we do not have enough gauge symmetries in order to gauge the 
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component e° of the standard vierbein to 1, whereas the component e° of the dual 
vierbein is constraint to be constant (that can be chosen as 1) by eq. (I2.2ip . Then 
the most general ansatz is as follows (where a, b, r, s, a, r, s, are functions of t): 

e a M := diag(6, a, a, a) 

ab Oi _ m ij _ ij 
U /j, ' j - r0 j> U k - Se k 

:= diag(l, a, a, a) 

~.ab . ~0i _ ~« ~ij _ ~Jj 

r(l) c. .(1) j r(l) * _ jW-. fc 

aft ■ ^ iO — 1 °i 5 ^ ij — 1 £ ij 

L^\ b :L { V = (4-1) 
(concerning the Lagrange multipliers, we have taken care of the constraints (I2.15P 
and 

Plugging these ansatze into the equations (and constraints) of section 2, we 
obtain from eq. (12.191) (where dots denote time derivatives): 

a = -r , (4.2) 

from (12.141) with (abc) = (iOj) and (ijk): 

a 2 r = —a 2 s , (4.3a) 

a 2 Sf = -a 2 r , (4.3b) 
from fl2T3H with (a6) = (iO) and (ij): 

= 45 d (£ {1) a 2 ) , (4.4a) 

2a d Q U 2 ) = . (4.4b) 



These equations are used to simplify some of the equations below, notably to 
eliminate £^ 2 \ 

For the stress energy tensor Tm we will make, to start with, the general ansatze 
Tmo > T Mi 3 = Tm,s , (4.5) 

and for T M 

Tm, ' = <V T M ,s ■ (4.6) 
Note that T M0 ° does not contribute to eq. ( 12. 29ft . As stated before, this does 
not constitute a paradox, since T M0 ° has to be determined by the conservation law 
that assumes, in the present context, the form 

a 
> 

a 
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Tmoo + 3? (T M00 - f MtS ) = . (4.7) 



Next eq. (12^5]) gives, for (ag) = (00) and (ij), 

3 (r 2 - s 2 ) = k a 2 T M oo , (4.8a) 



- 2r + - (r 2 - s 2 ) + Ak (V x) s - £ (1) f) = « a& T M ,s • (4.8b) 

Eq. (12.271) gives, for (abfi) = (Oij) and (ijk) (the other index combinations just give 
= 0): 



6 ab 

s = 4 9 (W) . (4.9b) 
Finally eq. (12.291) gives 

4a - - 29 «9 (£ (1) a 2 ) = -a 2 f MjS . (4.10) 

No further equations can be derived, and we are left with indeed 10 equations 
flOl-flOL fl4~8l-( |4TT0l) for 10 functions in the ansatz KTh . 

First it can be checked that, in the absence of dual matter {Tm,s = in eq. 
(I4.10p ). we can put s = £^ = I^ 1 = (and a = const.), and equations (14.81) collapse 
to (using (T4"IE1) ) 



a 2 



a 2 b 2 



k Tmoo , (4.11a) 



d a 2 ab m . 

2 ^ + ^- 2 ^ = kT ^- ( 4 - llb ) 
ab z a z b z ab 6 

These equations are invariant under time-like diffeomorphisms that allow the 
gauge b(t) = 1, after which they turn into the standard FRW equations that enforce 
the conservation of the standard stress energy tensor. 

Next we have analysed the system for arbitrary standard cosmological constant, 

Tmoo = T M>S = A , (4.12) 

and arbitrary dual cosmological constant, 

f M ,s = A • (4.13) 

The system of equations can be reduced by eliminating f using (14.21) . s using 
(14 .3b j) . and introducing 

4 1} = a 2 = const. (4.14) 
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that solves f!4.4bj) . Next, r can be eliminated using (I4.8al) . and b using f!4.9bj) . We 
are left with 4 equations for a, a, s and where the maximal time derivatives are 
d, d, s and i^: 




jW = - 2 4 i + % - ««? - - + iS (4.15) 

cr \ k a z / a gt da 6 I 



where one has to replace 

d 2k? ( (0 / .l ...n 
0= 1 z^— (aa — aa , (4.16a) 



r rcr 



r = ±^s 2 + ^a 2 A . (4.16b) 

This system can be brought into normal form (i.e. be solved for d, d, s and i^ 1 ** 
after the replacements (4.16)) which is suitable for analytic and numerical stability 
analyses. 

Remarkably we found for a wide range of initial conditions, and arbitrary 
cosmological constants A and A (and I* )> an asymptotically stable (constant) 
solution for the dual scale factor d: 

d — > d 
a — > ao/t 
s -> s 

-^t 2 A/4 (4.17) 

where the constants do, ao and so depend on the initial values. Another stable 
solution is given by 

d — > tao 

a — > ao/t 

s -> s 

-^t 2 A/24 (4.18) 
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that is separated from (I4.17P mainly through the initial values for a and Note 
that the A-dependent values of £^ are assumed dynamically. 

Provided that we identify our known matter with the dual matter of the model, 
these solutions come very close to a solution of the CCP, since the dual scale factor 
remains asymptotically constant (or increases just linearly in t) for arbitrary A, 
without fine tuned initial conditions. 

Let us compare the above solutions to the equations of the previous section 3, in 
particular to the points i) to v) just before eq. (3.1), in order to investigate at which 
point the impact of A on the dual scale factor a gets lost. First, the impact of A on 
L^\ b c in the form of £^ is evident. The crucial equation is eq. (4.9b), where the 
impact of £^ ~ t 2 on s (which corresponds to torsion lu t ) is cancelled through the 
decay of a(t) ~ l/t, up to subleading terms that allow for s — ► s . (Remarkably this 
does not imply that the "original" universe is contracting: from (4.16a) one obtains 
b(t) ~ bo/t 2 , hence t does not correspond to the cosmological time in this universe. 
Its cosmological time is rather given by t' = —b /t, in terms if which the signs of 
all components of e a M change and a(t') increases as \a(t')\ ~ ^t'. Note furthermore 
that increasing \t'\ corresponds to decreasing \t\, i.e. the relative arrows of time are 
reversed.) 

However, before this can be considered as a fully acceptable solution of the 
CCP, the following tasks have to be performed: 

a) dual matter has to be added in order to check, whether the known part of the 
evolution of our universe can be reproduced in the dual sector, 

b) a weak field expansion around such a solution has to be performed in order to 
see whether the (linearized) Einstein equations for the dual metric do not deviate 
too much from its standard form. 

We have performed a preliminary analysis in this direction by adding a 
cosmological perturbation AT to T = A. In the case of the validity of the standard 
Einstein equations this perturbation should induce a perturbation 

25 = kAT (4.19) 

of the dual scale factor a (cf. (4.11b)). The good news is that, neglecting terms 
of relative order ~ t -1 , the induced perturbation of a can indeed be written in the 
form (4.19). The bad news is that in both cases k is time dependent as k ~ t" 1 in 
the case of (4.17), even k ~ t~ 3 in the case of (4.18). (A similar problem persists 
for the "original" universe, which differs also from a de Sitter universe for arbitrary 
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A, unless A = and fine tuned initial conditions s = a = = are used, i.e. the 
CCP seems to be solved also here: however, the t 2 dependent value for plugged 
into the equations of motion (2.23), shows that the effective gravitational coupling 
for the "original" universe - its dependence on the Lorentz vector n d - is also time 
dependent.) Such time dependent gravitational couplings seem to be in conflict with 
the cosmological standard model, hence further studies of possible modifications of 
the model are required. 

5. Summary and outlook 

We have constructed an action including dual gravitational fields, without 
adding new degrees of freedom. An amazing feature of the model is that some kind 
of matter can be coupled to the original gravitational fields (and "see" a space-time 
geometry described by the original metric), whereas another kind of matter can be 
coupled to the dual metric and, hence, propagate on a generally different space-time 
geometry. In spite of the manifest breaking of Lorentz symmetry in the action, 
the reactions of the two different space-time geometries to the corresponding stress 
energy tensors coincide with Einstein's gravity in both cases, to lowest order in a 
weak field expansion around Minkowski space-time. 

Beyond lowest order, this phenomenon persists only for the standard gravi- 
tational fields (in the absence of dual matter): Whereas fluctuating standard 
gravitational fields imply fluctuating dual gravitational fields from eqs. (12.141) and 
f l2TT9|) . a vanishing f Ma b allows for vanishing LW and from f[2T29]) and fTOTD . 
as a consequence of which eqs. (12.231) and (12.271) turn into the standard equations 
for the standard vierbein and connection. 

The opposite statement is not true; an attempt to integrate out the standard 
gravitational fields (even for T Ma b = 0) will generate non-local effective interactions 
for the dual fields, that are moreover expected to break Lorentz symmetry manifestly 
in the form of (positive) powers of d^/n. Hence it is an open question up to now, 
whether our world could be identified with the dual sector of the model. 

This question is not purely academic, since the cosmological evolution of the 
dual sector differs dramatically from standard cosmology: A cosmological constant 
A does not imply an exponential increase (with cosmic time) of the dual scale factor 
a, but can be "absorbed" completely into a time dependence of the original scale 
factor. As stated at the end of section 4, such a solution of the cosmological constant 
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problem requires further investigations. 

Furthermore, various ways to generalize the model could be studied: 

a) As stated near the end of section 2, the fixed Lorentz vector n a could be 
replaced by a field n a {x) that breaks Lorentz symmetry spontaneously. The present 
approach to gravitational 5-duality would then be similar in spirit to the PST 
approach [23] (reviewed in [H]), that has been applied to d — 10/11 supergravities. 

b) In higher dimensional, e.g. 5 dimensional, space times the dual of R °* 
would be, after contracting the 5d epsilon tensor with the Lorentz indices, a tensor 
R abc ^ u that can be interpreted as a field strength of a field D abc ^ (antisymmetric in 
[a&c]). As before, this field could be fixed completely by an axial gauge condition and 
the duality constraint contracted with = n 5 . Subsequently the field strength of 
uj ab ,, can be obtained from n r R abc ,„,, and e a ,, from the torsion tensor contracted with 

[J, L ' flls I fj, 

n a as before (imposing the axial gauge condition on all these fields). A systematic 
procedure for arbitrarily high dimensional space-times could be developed along 
these lines, and it would be interesting to study the duals of brane universes in such 
models. (The 4d metric that is S-dual to 2-brane universes could be investigated 
already in the present model). 

c) A weak point of the present model is that, in spite of the introduction of 
dual gravitational fields, we did not manage to make additional symmetries (as of 
the Ehler's type [24]) manifest, even after dimensional reduction along a coordinate 
along n M . The technical problem here is that the standard t independent dimensional 
reduction ansatz for the metric does not coincide with our axial gauge conditions 
on e a ,, and oj ab ,,. 

It may then be advisable to carry out the essential steps of the concept presented 
here - the introduction of dual gravitational fields together with non-covariant gauge 
conditions and partial (non-covariant) duality constraints - in terms of different 
variables as non-linear realizations of gravity [25], or to give up the axial gauge 
condition on e 0, (which would complicate the analysis considerably, once n M can no 
longer be assumed to be constant). 

In view of the interesting properties of the present model we believe that these 
various open questions merit corresponding studies. 
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